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Abstract 

In the framework of the rigged Hilbert space, unstable quantum systems associated 
with first order poles of the analytically continued S-matrix can be described by Gamow 
vectors which are generalized vectors with exponential decay and a Breit-Wigner en- 
ergy distribution. This mathematical formalism can be generalized to quasistationary 
systems associated with higher order poles of the S-matrix, which leads to a set of 
Gamow vectors of higher order with a non-exponential time evolution. One can define 
a state operator from the set of higher order Gamow vectors which obeys the exponen- 
tial decay law. We shall discuss to what extend the requirement of an exponential time 
evolution determines the form of the state operator for a quasistationary microphysical 
system associated with a higher order pole of the S-matrix. 



1 Introduction 

Gamow vectors in rigged Hilbert spaces (RHS) were introduced [[L| to describe resonances, 
since they possess all the features usually attributed to resonance states, in particular, they 
obey an exponential decay law and have a Breit-Wigner energy distribution. They can be 
constructed from the first order poles of the analytically continued S-matrix on the second 
Riemann sheet of the complex energy plane, and are generalized eigenvectors of a self-adjoint 
Hamiltonian with complex eigenvalues (energy and lifetime) . Considering the S-matrix poles 

* Dedicated to L. C. Horwitz on the occasion of his 65 th birthday. 
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on the second Riemann sheet of higher order, e.g. order r, one can construct in the same way r 
Gamow vectors of orders k = 0, 1, 2, . . . , r— 1. These higher order Gamow vectors are Jordan 
vectors of degree k + 1, [Q] and they span an r-dimensional subspace of M. ZR C $ x , where $ x 
is the space of generalized vectors (functionals) of the Rigged Hilbert space $ C 7i C $ x . 
For the self-adjoint Hamiltonian H one obtains in the space $ x a matrix representation 
which, when restricted to the subspace M. ZR , form a Jordan block of degree r. || 

Jordan block matrices for the Hamiltonian were already investigated by Katznelson || 
who based his work on the well established ideas of Horwitz ||. However, their Hamil- 
tonians are non-self-adjoint, and their matrices are a generalization of the standard finite 
dimensional complex diagonalizable matrices for "effective" non-hermitian Hamiltonians ||. 
This approach cannot be implemented in a quantum theoretical framework using the Hilbert 
space. In contrast, the Gamow- Jordan vectors that are derived from higher order S-matrix 
poles have, like ordinary Gamow vectors representing Breit-Wigner resonances, a natural 
place in the rigged Hilbert space formulation of quantum mechanics. 

In section |2|, we will give a brief introduction to rigged Hilbert spaces, ordinary and higher 
order Gamow vectors, and their time evolution. In section 0, we will define the exponentially 
decaying higher order Gamow state operator which was conjectured in reference ]7) and 
discuss its generalization. We argue that these higher order Gamow state operators are of 
this form because (with certain qualifications) this is the only way to reconcile them with 
the exponential law. 

2 Resonances, Rigged Hilbert Spaces, 
and Gamow Vectors 

Conventionally, a resonance is described by a pair of first order poles of an analytically 
continued S-matrix jl]. These poles lie on the second sheet of the two-sheeted Riemann 
surface of complex energy at the complex conjugate positions z R = E R — iT/2 and z* R = 
Er + iT/2, where E R is the resonance energy and h/T is the lifetime of the resonance. The 
pole in the lower half-plane can be associated with a state that decays exponentially (for 
t > 0), while the pole in the upper half-plane can be associated with a state that grows 
exponentially (for t < 0). While in Hilbert space there are no such states that have an 
exponential time evolution and distinguish a direction of time, in the rigged Hilbert space 
(RHS) one describes such states by antilinear continuous functionals, the Gamow vectors. 
The rigged Hilbert space consists of a triplet of spaces ||: 

$CftC$ x . (2.1) 
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The best known example is the RHS where $ is the space (Schwartz space) of "well-behaved" 
functions, i.e. functions, that have derivatives, which are all continuous, smooth, and rapidly 
decreasing. TC is the space of Lebesgue square integrable functions, and <3> x (the space of 
continous antilinear functionals on $) is the space of tempered distributions. 

In quantum theory, if one distinguishes between preparations and registrations §|, one 
can further specify the RHS, and is led to a pair of RHS's: one for the preparations and one 
for the registrations []T0| . A scattering experiment can be subdivided into a preparation stage 
and a registration stage. The in-state + that evolves from the prepared in-state m outside 
the interaction region is determined by the preparation apparatus (the accelerator). The 
out-state ip~, detected as the "out-state" ip ont outside the interaction region, is determined 
by the registration apparatus (detector). According to the physical interpretation of the 
RHS formulation, "real" physical entities connected with the experimental apparatuses, e.g. 
the ensemble describing the preparation apparatus (the energy distribution of the 

beam) or the observable describing the registration apparatus (the energy resolution 

of the detector) are described by the "well-behaved" vectors 0, ip G $. One denotes the 
space of state vectors <p + by $_ and the space of "observable vectors" ip~ by $ + , where 
$ = $_ + $+ and $_ fl $ + 7^ 0. <£>_ is the space of "well-behaved" Hardy class vectors from 
below and $ + is the space of "well-behaved" Hardy class vectors from above, i.e. they fulfill 
even stronger properties connected with their analytic continuation, than the elements of 
the Schwartz space [[0]]. We will call these elements "very well-behaved" vectors. In place 



of the single rigged Hilbert space we therefore have a pair of rigged Hilbert spaces: 

(f) + G C 7i C for ensembles or prepared in-states, (2.2) 

G $+ C H C for observables or registered "out-states". (2.3) 

Here the Hilbert space 7i is the same for both triplets. 

On the level of $ or 7i one cannot talk of single microsystems, and there are no math- 
ematical objects in Hilbert space quantum mechanics to describe a single microsystem or a 
single experiment which prepares and observes a single microsystem. Still, it is intuitively 
attractive to imagine that the effect by which the preparation apparatus acts on the regis- 
tration apparatus is carried out by single physical entities, the microphysical systems. The 
energy distribution for a microphysical system does not have to be a "well-behaved" func- 
tion of the physical values of energy E. Hence, for the hypothetical entities connected with 
microphysical systems, like Dirac's "scattering states" |p) or Gamow's "decaying states" 
\E - iT/2), the RHS formulation uses elements of $ x , $*, and $ x . 

The decaying Gamow vector associated with the complex energy zr = Er — iT/2 is 
a continuous antilinear functional over the vectors ip~ G $+, and its time evolution is 
defined for times t > 0. They are generalized eigenvectors of (the extension of) a self-adjoint 
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Hamiltonian H with complex eigenvalue zr 



H x \z R ) = z R \z R ); z R = E R -iY/2 (2.4) 

(where the x denotes the conjugate operator acting on the functionals). These Gamow 
vectors can be obtained from the first order poles of the analytically continued S-matrix at 
the position zr on the lower half-plane of the second Riemann sheet, 

S(z) = h analytic terms . (2.5) 

z — z R 



Here a_i can determined from the unitarity of the S-matrix to be — iT fI2 |. In the same way 
one can define the exponentially growing Gamow vector associated with the complex energy 
z* R = Er -Mr/2 , which is a continuous antilinear functionals over the vectors <p + G 

The n th order Gamow vector associated with the complex energy zr is a generalization 
of the Gamow vector above in the following way: 

1) There are r generalized vectors (functionals) of order n = 0, 1, . . . , r — 1. They are 
associated with a pole of the order r at the position z R on the second Riemann sheet of the 
analytically continued S-matrix, 

S(z) = h - — — + • • • + - ^— — h analytic terms (2.6) 

z - z R (z- z R y (z - Zr) t 

where the coefficients a_ n _i, n = 0, 1,..., r — 1 can be determined from the unitarity of 
the S-matrix and are complex numbers with dimension [energy] n+1 . 

One starts from the S-matrix elements, i.e. the matrix elements of the incoming prepared 
state m G $_ at t — > — oo and the outgoing detected observable ijj out G at t — >• oo 

(^ out ,50 in ) = (7/; out ,fi~ t fi + in ) = (Q-ij out ,Q + (P in ) 

POO 

= (vr,0 + )=/ dE(^-\E-)S(E)( + E\ ( f> + ) , (2.7) 
Jo 

where fi ± are the M0ller wave operators. One can now deform the contour of integration 
from the cut (positive real axis and spectrum of H), into the lower half-plane of the second 
sheet. Then one obtains a "background integral" term, independent of the poles, along the 
negative real axis of the second sheet and a residue term for the pole: 



0/r 



/ dE(^-\E-)S(E)( + E\(f) + ) + residue at z R . (2. 

J -OOrr 
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For every residue term of the S-matrix pole at z R one obtains: 
residue at z R = nT^dz^^^^^^ 

n=0 

-27rm_ ra _i ^ / n\ d n ~ k d k + , 

= e— — M*Js^ 1 } ^ ( " 10 } 

n=0 fc=0 V 7 

r-1 



(2.9) 



n=0 ' fc=0 ^ 7 



where (V'~l- 2; ij) < ' n_A: ' ) is t ne ( n ~~ fc)-th derivative of the "well-behaved" (continuous, analytic, 
smooth, rapidly decreasing) function {ip~\z~) G 7^ 2 PI S at the position ;zr, and is therefore 
analytic in the lower half-plane, and ^ k \ + z R \(j) + ) is the fc-th derivative of the "well-behaved" 
analytic function ( + z\(p + ) G 1-L 2 _ D S at the position z R in the lower half-plane. The S- 
matrix ( [2.6|) is thus associated with a set of r generalized vectors, 

\z R f\ \z R f\..., \z~ R f-V (2.10) 

where \z R f- ^ is the ordinary Gamow vector. (The same generalization can be carried out for 
the growing higher order Gamow vectors with complex energy z* R which we shall not discuss 
here.) 

2) It can be shown that the higher order Gamow vectors are generalized eigenvectors of 
the self-adjoint Hamiltonian H in the following sense 

{Hi,-\z R ) = &-\H*\z R p = z R ^-\z R p 

m~\z R p = (r\H*\z R p = z R {r\z R p + (r\z R p 

'■■ (2.11) 

(ifx/r|z-)« = (ip-\H x \z R p = z R (ij-\z R f k) + k{ilj-\z R f- 1) 

(Hr\z R f- l) = &-\H*\z R f-V = z R (^-\z R j r -V + (r-l)(^~\z R f-V 
Omitting the arbitrary ip~ one writes this as 

H x \z R f k) = z R \z R )W + k\z R f k - 1 ) . (2.12) 
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From this it follows that H x is a Jordan operator of degree r, and the k-th order Gamow 
vector \Zfif k > is a Jordan vector of degree k + 1 ||, i.e., 

(H x - z R ) {k+1) \z R fV =0 and (F x - z^^l^f^ ^ . (2.13) 

3) The time evolution of the decaying fc-th order Gamow vector is given by [2] 



e~* HXt \z 



J*> = e+** (t) (-*t) k - p \z R r\ t > 0. (2.14) 
p=0 ^ ' 



The fc-th order Gamow vector evolves into a superposition containing Gamow vectors of the 
same and all lower orders and the space spanned by the set (2.10) is thus invariant under 
the action of the time evolution operator. It follows from (2.14) that the time evolution 
of the dyadic product of a k-th order Gamow vector with an m-th order Gamow vector, 
\ z R~f kS> ^ m \~~ z r\i has terms with additional powers of t multiplying the overall exponential 
factor e~ rt . Such non-exponential time evolution is really no surprise in view of the previous 
results [14]. A more surprising fact is that certain linear combinations of the dyadic products, 
e.g., 

W (n) = — Z)U) l*H> (fc) (n ~*Wl> for fixed {0,1,- ■• ,r-l}, (2.15) 
n ' fc=0 ^ ' 

have a purely exponential time evolution [7]: 

W in \t) = e- iHXt W {n) (0)e iHt = e- rt W (n \0), for t > 0. (2.16) 



3 Form of the Exponentially Decaying Operators 

We now wish to determine to what extent the requirement of exponential time evolution 
restricts the form of an operator constructed as a linear combination of dyadic products of 
vectors in M. ZR - The most general linear combination of dyadic products of vectors in M. ZR 
is given by 

r— 1 r— 1 
k=0 m=0 
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with arbitrary coefficients B m k . We will show that W decays according to the pure expo- 
nential e~ Tt if and only if the coefficients are restricted by 

Bmk={ y^B^o, for* + m<r-l, (32) 

0, for k + m > r — 1, 

where the coefficients B k+m0 remain arbitrary, that is, if and only if W is restricted to the 
form 



n=0 k=0 



^ = £ 5 n,o£(? ) \znf k) {n - k) (-*R\ (3-3) 



with arbitrary coefficients B nfi . Since the vectors |%)^ have the dimension [energy]" 5 ~ fc , the 
sums over in (3.3) have the dimensions [energy]" 1 "™ so that, if 2irTW is to be dimensionless, 
then the coefficients B n $ must have the dimensions [energy]" as in (2.15). 

For the proof it is convenient to consider the most general linear combination of dyadic 
products \z~ R )^ ( m \~z R \ for which the sum n = k + m of the orders does not exceed a finite 
integer j; this is given by 



3 n 



n=0 k=0 



with arbitrary coefficients A njk . The operator (3.1) is obtained as a special case of the 
operator (3.4) by setting 

, , . / 0, for n - k > r - 1, , , 

J = 2(r-1), ^ = ^-M=( 0)for fc>r _i. (3-5) 

The time dependence of W^j) is given, using (2.14), by 

W U) (t) = (e^r%(0)e !ffi = tE^(e !ffi ) x |%f )M) (-^|e !ffi 

n=0 k=0 

= «- n EEEE ^ (?) (""*) («)-*- 1^> (,) ( -T*«i. 

ri=0 fc=0 2=0 m=0 V / V / 

Changing the order of the summations, 

j n k n—k j n n n—k j j n n—k j j n—l n—m j j—l j n—m 

EEEE=EEEE=EEEE=EEEE=EE E E- 

n=0 k=0 1=0 m=0 n=0 1=0 k=l m=0 1=0 n=l k=l m=0 1=0 n=l m=0 k=l 1=0 m=0 n=l+m k=l 
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allows the dyadic products, which are linearly independent operators, to be factored out of 
the sums over terms in which they appear as common factors: 

w m (t) = »- r 'EE E E M , ) (" ) (-«)*-'(«)"-'-"• < m >r^i 

/— m—rin — lJ-m h — 1 \ / \ / 



£=0 m=0n=l+m k= 
j j-l j n-m 



'k\ ( n — k 
m 



I ••\k— I / •j.sn—k—m 
(-It) (it) 



J J~ 1 J lb — III / / \ / 

e ""EEi^> ( ' , "" ) <" 2 *i E Emj( 

1=0 m=0 n=l+m k=l V / V 

j j~ l j n-m /U\ / U\ 

^ r 'EE \*-«)<" '-T'-i E (*)"-"-' E M?) (T„ ) <-D 



The operator W(j)(t) will decay according to the pure exponential e if and only if all 
terms involving additional powers of t cancel. All terms involving additional powers of t will 
cancel if and only if the coefficients A n k satisfy the conditions 



k\ (n — k 



( I G {0, • • • ,j-l}, 

= Vij;](" ")(-l) k - 1 for { mG{0 ,... ?J -l-/}, (3.6) 
k=i \ L J\ m J { ne { m + i + i r .. j}, 

The simplest of these conditions are those for which n = m + I + 1, i.e., those for which 
m — n — I — 1, because they are the only conditions that involve sums over only two values 
of k: 



k=i 



Z J f ° r i JG{0,...,n-l} 



or, replacing Z with fc — 1, 



n " fc + 1> \ -4 f M for ( ne{l,--- ,j}, 

n-k ) - An Ak-l tOT 1 fcG{l,...,n} 



or, using the definition (™) = fc!(ra "i fc)! , 



_ ( w -fc + l)!(fc-l)! f nelV-j}, 

A "' fc " (^kjkl An ' k ~ l tOT 1 fc G {!,•••, n}. 



8 



These conditions relate pairs of coefficients A n ^ having the same values of n and successive 
values of k. For fixed n G {1, 2, • • • , j}, they may be used recursively to show that A Ujk must 
equal A„ i0 multiplied by (™) : 



A 



n,k 



(n-k+l)\(k-l)\ 



(n-k)\k\ 



{n-k+2)\(k-2)\ 
{n-k + l)\{k-l)\ 



n!0! 



{n-k)\k\ 



A n ,o 



A n ,o 



for 



(w-l)!l! 
(n-2)!2! 

n e {!,-•• 

fc G {!,-•• ,n}. 



n!0! 



(n-l)!l! 



Substituting this result into the full set of conditions (3.6), using the identity 



n — k 
m 



n\ I n — m\ f n—m — l 



rn 



I 



k-l 



and then using the binomial formula gives 




A 



l n,0 



k=l 

n 
rn 



n — k \ 

'_1 f 1 

m ' 



A 



71,0 



A 



n,0 



= A 



71,0 



n—m 
I 

n—m 
I 

n—m 
I 

n—m 
I 



E 

k=i 

n—m— l 

E 

fc-Z=0 

(i-ir 

0' 



n—m — l 
k-l 

n — m — l 
k-l 



-1) 



fe-j 



-^n—m—k / ! 



for 



Ze{0,...,j-l}, 

nG {m + Z + !,-•• ,j}, 



n,0 



(3.7) 



which shows that the remaining conditions are automatically satisfied by (3.7) without plac- 
ing any further conditions on the coefficients A n>0 . The coefficients A n> o, for n G {1, • • • ,j}, 
and also the coefficient Aq^, remain completely arbitrary. 

We conclude that a linear combination of dyadic products \z^f k ^ ^ m \~ZR\ for which the 
sum n = k + m of the orders does not exceed j decays according to the pure exponential 
e~ rt if and only if it is of the form 



3 n / \ 

E^.oE(>i> l ' )< "-' , w 

ra=0 k=0 ^ ' 



(3.8) 



9 



with arbitrary coefficients A n . Equation (3.3), i.e., (3.1) with the restrictions (3.2), follows 
by applying the special case (3.5) to the result (3.8). This concludes the proof. 



4 Summary and Conclusion 

In the conventional Hilbert space quantum mechanics, resonance states cannot be described 
by state vectors. Therefore the most common definition of a resonance is as a pole of the 
S-matrix at the complex energy zr = Er — iT/2 (actually, a pair of poles at Er ± iT/2). 
There is no theoretical reason to exclude poles of order higher than one, and poles of any 
order will lead to the typical resonance phenomena for cross-section and phase shift, in 
particular, the time delay and therewith formation of a quasi-stationary state |T2|. But 



higher order poles have been scorned, because they were somehow associated with a time 
evolution that in addition to the exponential had also a strong (of order h/T) polynomial 



time dependence |14j] for which there exists no experimental evidence. In the rigged Hilbert 
space formulation of quantum mechanics, a vector description of resonances is possible, and 
Gamow vectors and Gamow states were defined from the first order pole of the S-matrix 
element. Their Hamiltonian time evolution was shown to be exactly exponential and given 
by a semigroup IT3"] . This procedure was generalized to an r-th order pole at z = zr 
which led to the definition of r Gamow vectors of order k — 0, 1,..., r — 1. The 
Gamow vectors of order k > 1 were shown to have indeed a polynomial time evolution in 
addition to the exponential J|. However, one can find a non-reducible state operator in 
the r-dimensional space M. ZR spanned by the r Gamow vectors, e.g. the operator (|2.15|), 



which has a purely exponential time evolution [7J. In the paper at hand we investigated 
the question, to what extend the requirement of a purely exponential time evolution with 
a lifetime r = h/T = h/lmzR determines the form of the state opeerator. From rather 
plausible assumptions we concluded that this must be a "mixture" of operators like ( [2.15 ) 
with arbitrary coefficients A n $ as given by 
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